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14.  ABSTRACT 

For  high  Mach  number  flows,  M>  4,  piston  theory  has  been  used  to  calculate  the  pressures  on  the  surfaces  of  a  vehicle.  In  a  two- 
dimensional  flow,  a  perpendicular  column  of  fluid  stays  intact  as  it  passes  over  a  solid  surface.  Thus,  the  pressure  at  the  surface  can  be 
calculated  assuming  the  surface  were  a  piston  moving  into  a  column  of  fluid.  In  this  work,  first-order  piston  theory  is  used  to  calculate 
the  forces,  moments,  and  stability  derivatives  for  longitudinal  motion  of  a  hypersonic  vehicle.  Piston  theory  predicts  a  relationship 
between  the  local  pressure  on  a  surface  and  the  normal  component  of  fluid  velocity  produced  by  the  surface's  motion.  The  advantage 
of  piston  theory  over  other  techniques,  such  as  Prandtl-Meyer  flow  or  Newtonian  impact  theory,  is  that  unsteady  aerodynamic  effects 
can  be  included  in  the  model.  The  unsteady  effects,  considered  in  this  work,  include  perturbations  in  the  linear  velocities  and  angular 
rate.  This  provides  a  more  accurate  model  that  agrees  more  closely  with  models  derived  using  computational  fluid  dynamics  or  those 
derived  by  solving  Euler  equations.  Additionally,  piston  theory  yields  an  analytical  model  for  the  longitudinal  motion  of  the  vehicle, 
thus  allowing  design  trade  studies  to  be  performed  while  still  providing  insight  into  the  physics  of  the  problem. 
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I.  Abstract 

For  high  Mach  number  flows,  M  >  4,  piston  theory  has  been  used  to  calculate  the  pressures  on  the 
surfaces  of  a  vehicle.  In  a  two-dimensional  inviscid  flow,  a  perpendicular  column  of  fluid  stays  intact  as  it 
passes  over  a  solid  surface.  Thus,  the  pressure  at  the  surface  can  be  calculated  assuming  the  surface  were  a 
piston  moving  into  a  column  of  fluid.  In  this  work,  first-order  piston  theory  is  used  to  calculate  the  forces, 
moments,  and  stability  derivatives  for  longitudinal  motion  of  a  hypersonic  vehicle.  Piston  theory  predicts  a 
relationship  between  the  local  pressure  on  a  surface  and  the  normal  component  of  fluid  velocity  produced 
by  the  surface’s  motion.  The  advantage  of  piston  theory  over  other  techniques,  such  as  Prandtl-Meyer  flow, 
oblique  shock,  or  Newtonian  impact  theory,  is  that  unsteady  aerodynamic  effects  can  be  included  in  the 
model.  The  unsteady  effects,  considered  in  this  work,  include  perturbations  in  the  linear  velocities  and 
angular  rates,  due  to  rigid  body  motion.  This  provides  a  more  accurate  model  that  agrees  more  closely  with 
models  derived  using  computational  fluid  dynamics  or  those  derived  by  solving  Euler  equations.  Additionally, 
piston  theory  yields  an  analytical  model  for  the  longitudinal  motion  of  the  vehicle,  thus  allowing  design  trade 
studies  to  be  performed  while  still  providing  insight  into  the  physics  of  the  problem. 

II.  Introduction 

In  the  I980’s,  the  National  Aerospace  Plane  (NASP)  program  commenced,  with  its  goal  being  a  feasibility 
study  for  a  single-stage  to  orbit  (SSTO)  vehicle,  which  was  reusable  and  could  take  off  and  land  horizontally. 
The  NASP  was  to  be  powered  by  a  supersonic  combustion  ramjet  (scramjet)  engine.  Although  this  program 
was  cancelled  in  the  I990’s,  a  great  deal  of  knowledge  was  gained  and  it  spawned  future  programs,  including 
the  hypersonic  systems  technology  program  (HySTP),  initiated  in  late  1994,  and  the  NASA  X-43A.  The 
HySTP  program’s  goal  was  to  transfer  the  accomplishments  of  the  NASP  program  to  a  technology  demon¬ 
stration  program.  This  program  was  cancelled  in  early  1995.  The  NASA  X-43A  set  new  world  speed  records 
in  2004,  reaching  Mach  6.8  and  Mach  9.6  on  two  separate  occasions  with  a  scramjet  engine.  These  flights 
were  the  culmination  of  NASA’s  Hyper-X  program,  with  the  objective  being  to  explore  alternatives  to  rocket 
power  for  space  access  vehicles. 

With  renewed  interest  in  space  operations  worldwide,  there  is  a  renewed  interest  in  hypersonic  aerody¬ 
namics  research.  The  scramjet  engine  will  likely  play  a  major  role  in  future  hypersonic  vehicles.  Ulinke  a 
conventional  turbojet  engine,  a  scramjet  engine  does  not  use  high  speed  turbomachinery  to  compress  the  air 
before  it  reaches  the  combustor.  Instead,  it  relies  upon  the  rise  in  pressure  across  oblique  shock  waves  located 
in  front  of  the  inlet.  Furthermore,  the  flow  through  the  entire  engine  is  supersonic  in  contract  to  a  ramjet 
where  the  flow  speeds  are  subsonic  through  the  combustor.  On  configurations  like  the  NASP  and  X-43A, 
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the  underside  of  the  airframe  must  function  as  the  air  inlet  mechanism  and  the  exhaust  nozzle.  Therefore, 
integration  of  the  airframe  and  engine  are  critical  to  success  of  a  scramjet  powered  vehicle. 

Scramjets  could  be  used  as  part  of  a  multi-stage  launch  vehicle  that  would  include  multiple  propulsion 
systems  to  perform  a  mission.  The  factor  driving  research  towards  scramjets  and  away  from  rockets  is  cost; 
scramjets  would  substantially  lower  costs  because  it  is  an  airbreathing  engine.  Airbreathing  engines  don’t 
require  oxidizer  to  be  carried  by  the  vehicle,  hence  increasing  the  payload  and  reducing  the  quantity  of  fuel 
carried. 

Unsteady  aerodynamics  is  another  key  technology  in  the  development  and  optimization  of  future  hyper¬ 
sonic  vehicles.  The  combined  effects  of  a  slender  flexible  vehicle  travelling  at  high  speeds  and  subjected  to 
large  forces  may  lead  to  significant  unsteady  aerodynamic  effects.  Hence,  understanding  the  concepts  and 
consequences  of  time-dependent  aerodynamic  flows  is  critical  to  the  success  of  this  type  of  vehicle. 

Piston  theory  is  a  technique  which  has  been  used  for  years  to  determine  the  pressure  distributions  on 
an  air  foil/ vehicle,  when  the  Mach  number  is  sufficiently  high.  Lighthill^  discussed  the  application  of  piston 
theory  on  oscillating  airfoils  some  50  years  ago.  Ashley  and  Zartarian^  discuss  piston  theory  while  providing 
a  number  of  examples  of  the  application  of  piston  theory  to  specific  problems.  More  recently,  Tarpley^ 
discussed  the  computation  of  stability  derivatives  for  a  caret- wing  waverider  using  piston  theory.  Estimation 
of  stability  derivatives  also  requires  the  analysis  of  unsteady  flow  over  the  vehicle.^  Piston  theory  allows  the 
inclusion  of  unsteady  aerodynamic  effects  in  the  model  and  a  closed  form  solution  can  be  found  for  these 
unsteady  effects. 

In  this  work,  piston  theory  is  applied  to  a  hypothetical  2-dimensional  hypersonic  vehicle  powered  with 
a  scramjet.  Section  III  describes  the  vehicle  analyzed  in  this  work,  while  the  steady  forces  on  the  vehicle’s 
surfaces  are  calculated  in  Section  IV.  The  afterbody  effects  are  included  in  Section  V,  where  the  pressure 
distribution  and  force  due  to  exhaust  plume  are  evaluated.  The  control  surface,  which  is  a  single  elevator, 
is  examined  in  Section  VI,  analysis  of  the  flow  regions  is  performed  in  Section  VII,  total  body  forces  and 
moments  are  derived  in  Section  VIII,  while  stability  derivatives  are  calculated  in  Section  IX.  The  engine 
model  is  developed  in  Section  X  and  some  results  from  an  open-loop  simulation  are  provided  in  Section  XL 


III.  HSV  Model 

Figure  1  shows  the  hypersonic  vehicle  considered  in  this  work.^  The  vehicle  consists  of  4  surfaces:  an 
upper  surface  (surfaces  defined  by  points  cf)  and  three  lower  surfaces  (surfaces  defined  by  points  cd,  gh,  and 
ef).  All  pertinent  lengths  and  dimensions  are  in  units  of  feet  and  degrees,  respectively.  The  total  length  of 
the  vehicle  is  L  =  100/t  and  the  notation  for  lengths  is  Lf  =  length  of  the  forebody,  Lg  =  length  of  the 
engine  nacelle.  La  =  length  of  the  aftbody,  Lg  is  the  length  of  the  elevator,  Xf  is  the  distance  from  the  C.G. 
to  the  front  of  the  vehicle,  Xa  is  the  distance  from  the  C.G.  to  the  rear  of  the  vehicle,  Xcs  and  Zgg  are  the 
distances  from  the  C.G.  to  the  midpoint  of  the  elevator  in  the  x  and  z  directions,  respectively,  and  hi  is  the 
engine  height.  The  vehicle  lengths  are 

L  =  100  ft 
Lf  =  47  ft 
La  =  33  ft 
Ln  =  20  ft 
Le  =  17  ft 

Xf  =  55ft 

Xa  =  45  ft 
Xcs  =  30  ft 
Zcs  =  3.5  ft 
hi  =  3.5  ft 
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The  vehicle  angles  are 


^i,u  =  3° 

n.L  =  6°  (2) 

T2  =  14.41° 

Additionally,  the  vehicle  mass  and  moment  of  inertia  are 

Mass  =  300  slug  .  x 

Jyy  =  500, 000  slug  — 

and  the  mean  aerodynamic  chord  (c)  and  planform  area  {S)  are  defined  as 


c  =  L 

S  =  L^ 


(4) 


The  goal  is  to  apply  piston  theory  to  this  vehicle  to  determine  the  pressure  distribution  on  the  surfaces 
of  the  vehicle,  which,  in  turn,  can  be  used  to  evaluate  the  forces  and  moments.  The  pressure  on  the  face  of 
a  piston  moving  into  a  column  of  perfect  gas  is^ 


7-1  K 
2  clqo 


27 

7-1 


(5) 


where  the  subscript  ”  oo”  refers  to  the  steady  flow  conditions  past  the  surface,  Vn  is  the  velocity  of  the  surface 
normal  to  the  steady  flow,  ttoo  is  the  freestream  speed  of  sound,  and  P  is  the  surface  pressure.  Taking  the 
binomial  expansion  of  Eq.  5  produces 

T  =  1  + =  i  +  (6) 

PqO  7  1  ^  ^oo  ^oo 

Multiplying  through  by  Poo  and  using  the  perfect  gas  law  (P  =  pRT)  and  the  definition  of  the  speed  of 
sound  {o?  =  jRT)  yields  the  basic  result  from  first-order  linear  piston  theory 


P  —  -Poo  T  Poo^^oo^n 


(7) 


The  infinitesimal  force  due  to  the  pressure  is 


dF  =  -PdAn  (8) 

where  dA  is  a  surface  element  and  n  is  the  outward  pointing  normal.  Substituting  Eq.  7  into  Eq.  8  yields 

dF  =  (-Poo  -  PoodooVn)  dAn  (9) 


The  normal  velocity  can  be  computed  by  taking  the  dot  product  of  the  flow  velocity  over  a  surface  and  the 
outward  pointing  normal  for  that  surface.  Hence,  Eq.  9  becomes 

dF  =  (-Poo  -  PooCioo  [V  •  n])  dAn  (10) 


Equation  10  is  the  basic  result  upon  which  this  work  is  based.  Erom  this  equation,  it  is  seen  that  in  order 
to  compute  the  forces  acting  on  a  surface,  one  must  determine  the  properties  of  the  flow  past  the  surface 
(properties  behind  a  shock,  expansion  fan,  or  freestream),  the  velocity  of  the  surface  relative  to  the  airstream, 
V,  the  outward  pointing  surface  normal,  n,  and  the  surface  element,  dA.  The  work  that  follows  will  develop 
these  necessary  quantities  for  the  upper  surface  and  the  lower  surfaces  defined  by  points  (cd)  and  (gh). 
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Figure  2.  Hypersonic  Vehicle,  Oblique  Shocks,  and  Pressure  Regions. 


IV.  Vehicle  Surface  Pressure  Distributions  and  Forces 

To  compute  the  forces,  moments,  and  stability  derivatives,  consider  small  perturbations,  from  a  steady 
flight  condition  at  Moo,  in  the  velocities  u  and  w  and  the  rate  q.  On  the  upper  surface,  the  surface  is 
modelled  as  a  piston  moving  into  a  column  of  fluid  that  has  the  properties  of  the  flow  behind  an  oblique 
shock  wave  (an  oblique  shock  forms  if  a  <  Likewise,  on  the  lower  surface,  the  surface  is  modelled  as  a 

piston  moving  into  a  column  of  fluid  that  has  the  properties  of  the  fluid  behind  the  oblique  shock.  Figure  2 
shows  the  regions  of  interest. 

To  begin  the  development,  first  consider  the  upper  surface.  The  velocity  of  a  point  on  the  upper  surface 
due  to  the  velocity  and  rate  perturbations  is 

Vc/  =  ( Vi  cos  Ti^t/  +  1^)  i  +  ( Vi  sin  +  re)  ^  +  u;  x  rc/  (11) 

where  i,  k  are  unit  vectors  in  the  x  and  z  body  axes,  respectively,  u?  is  the  angular  rate  vector,  a  is  the  angle 
of  attack,  and  Vi  is  the  velocity  of  the  flow  behind  the  oblique  shock  wave  (see  region  1  in  Fig.  2).  For 
longitudinal  motion  only,  u?  =  qj  where  j  is  a  unit  vector  in  the  y  body  axis  direction.  In  Eq.  11,  is  the 
position  vector  of  a  point  on  the  upper  surface  given  by 

I*c/  =  TcfJ  +  +rcf,k  =  xi  +  tanrpc/  {x  -  Xf)  k 
—Xa  <  X  <  Xf 

According  to  Eq.  8,  a  normal  vector  to  the  upper  surface  is  also  needed.  The  upper  surface  outward  pointing 
normal  vector  is 

Ucf  =  sinri^t/i  —  cosri^uk  (13) 

Eor  the  lower  surface  defined  by  the  points  c  and  d  in  Eigure  1,  we  use  the  velocity  of  the  flow  after  the 
oblique  shock  to  obtain 

Vcd  =  (E2  costal  +  (-k2sinrqL  w)  k  ^  lv  x  Vcd  (14) 
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while  for  the  surface  defined  by  points  g  and  h 

^gh  =  (Vs  ^u)i^wk  +  X  Ygh 

where  Ycd  and  Ygh  are  position  vectors  of  a  point  on  the  lower  surface  given  by 

I*cd  =  TcdJ  +  Tcd.k  =  xi  -  tanrpL  {x  -  Xf)k 
Xf  —  Lf  <  X  <  Xf 


(15) 


(16) 


Ygh  =  TghJ  +  rgh^k  =  xi-\-  (L/tanryL  +  hi)k 
(xf  —  Lf)  —  Lji  <  X  <  Xf  —  Lf 
The  normal  vectors  for  the  lower  surfaces  are 

Ucd  =  sinrpL^  +  COSTyL^^^ 

Ugh  =  Ik 

Performing  the  cross  products  required  by  Eqs.  11,  14,  and  15  gives 

(jj  X  Ycf  =  gtanri^f/  {x  —  Xf)i  —  qxk 

u?  X  Ycd  =  — ^tanri^L  {x  —  x f)  i  —  qxk 
(jj  X  Ygh  =  g(I//tanTi,L  hi)i  -  qxk 
According  to  Eq.  10,  the  pressures  on  the  surfaces  of  interest  are 


(17) 

(18) 

(19) 

(20) 
(21) 


Pcf  =  Pi  Pldl  (Vc/  •  Ylcf) 

Pcd  =  P2  P  P2CI2  {Ycd  '  YVcd)  (22) 

Pgh  —  P3  T  PS^l^S  (Ygh  '  Ylgh) 

Substituting  the  results  of  Eq.  22  into  Eq.  8  gives 

dF cf  =  I  P I  picii  (Y cf  '  ^c/)}  dA-cfYVcf 

dF cd  —  {  P2  (Ycd  '  I^cd)}  dAcd^i-cd  (23) 

dF gh  =  {  P3  (Ygh  '  Ylgh)}  dAghYlgh 

Using  Eqs.  11,  14,  and  15  and  the  appropriate  normal  vectors  (Eqs.  13  and  18),  the  dot  products  in  Eq.  23 
become 

Vc/  •  Yicf  =  [u  P  gtanri^f/  {x  —  Xf)]  sinri^t/  —  [w  —  qx]  cosri^u 

Vcd  •  Ucd  =  [u-  gtanrpL  (x  -  Xf)]  sinr^L  P  [w  -  qx]  costal  (24) 

Vgh  'ngh  =  w  -qx 

Note  that  the  steady  terms  cancel  as  a  result  of  taking  the  dot  product.  Using  Eq.  24  in  Eq.  23  yields 


dFcf  =  (—Pi  —  PiCii  +  gtanri^f/  (x  —  Xf)]  sinri^t/  —  [w  —  qx]  cosri^t/})  dAcfUcf 
dFcd  =  (-P2  -  p2Ci2  {[u  -  gtanrpL  (x  -  Xf)]  sin  ryi,  P  [w  -  qx]  costal})  dAcd^icd  (25) 

dFgh  =  (-P3  -  p3Ci3  -  qx})  dAghUgh 

The  next  step  is  to  determine  the  upper  and  lower  surface  elements.  Note  that  the  vehicle  model  is  2- 
dimensional  with  unit  depth  into  the  page.  Hence,  the  upper  surface  element,  dAu  can  be  written  as 

dAcf=dLcf(l)  (26) 

where  dLcf  defines  a  length  of  interest  on  the  upper  surface  and  the  multiplying  factor  of  1  is  due  to  the 
vehicle’s  unit  depth.  The  surface  element  can  be  written  as 

dAcf  =  Y dx‘^  P  dz‘^(l)  (27) 
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From  Eq.  12, 


z  =  tanri^f/  {x  —  Xf)  dz  =  tanri^t/c^x 


Using  Eq.  28  in  Eq.  27  yields 

dAcf  =  dx‘^  +  tan^  Ti^udx‘^{l)  =  dx^J  1  +  tan^  =  dx  sec  ri^t/(l)  =  secri^t/dx 

Similarly,  the  surface  elements  for  the  lower  surfaces  become 


dAcd  =  secTi^Ldx 
dAgh  =  dx 

Now,  the  incremental  force  on  the  upper  surface  (see  first  entry  in  Eq.  25)  becomes 

dFcf  =  (—^1  —  pidi  {['^  +  ^tauTi^f/  {x  —  Xf)]  sinri^t/  —  [w  —  qx]  cosri^t/})  Ucf  secri^t/dx 
Using  similar  analysis,  the  incremental  forces  on  the  lower  surfaces  become 


(28) 

(29) 

(30) 

(31) 

(32) 


dFcd  =  {-P2  -  P2a2  {[u-  gtanri^L  (x  -  Xf)]  sinri^L  +  [w  -  qx]  costi^l})  UcdsecTi^Ldx  (33) 

dFgh  =  {-P3  -  Psaa  {w  -  qx})  Ughdx  (34) 

It  will  be  useful  to  separate  the  steady  forces  from  the  unsteady  forces  at  this  point.  In  Section  IX,  we 
will  model  the  unsteady  effects  using  a  stability  derivative  approach.  To  compute  the  steady  forces,  the 
components  of  the  incremental  forces  related  to  the  steady  flow  are  integrated  over  the  surface  of  the  vehicle. 
For  the  upper  surface,  this  becomes 


Fc/  = 


r^f 


{dFcj) 


r^f 


steady 


-Pi 


siuTi  [/i  —  cosTi  uk 


secTi  udx 


(35) 


while,  for  the  lower  surfaces,  the  forces  are 


f^f 

steady 


Fed  = 


r^f 


-P2 


rXf-Lf  POL 

'^gh  =  /  steady  ^  / 

J{xf-Lf)-Ln  J(xf-Lf)-Ln 

Eor  the  consideration  of  aeroelastic  effects  in  future  work,  it  will  prove  useful  to  split  the  upper  surface  force 

into  forebody  and  afterbody  components.  Thus,  Eq.  35  becomes 

ro 


(dPgh)  steady 


sinri^L^  +  costi^l^ 

Xf-Lf 


secTi^Ldx 
—Pskdx 


(36) 

(37) 


Fc/  =  /  -Pi 

J  —Xn 


siuTi^t/t  —  cosTi^t/^  secTi^t/dx  +  /  — 

-I  ’  do 


Pi 


siuTi^t/t  —  cosTi^uk  secTi^udx 


(38) 


'cfa 


'eff 


Performing  the  integrations  yields 


Fc/„  =  -PlXaSeCTi^u  siUTyt/t  -  COSTyf/^ 

Feff  =  —PlXf  SeCTi^u  siUTi^t/t  —  COSTi^uk 


~  ^efa  T  ^efa  k 

—  ^eff'^  T  ^eff  k 


(39) 


where  Xcf^ ,  Zef^  are  the  components  of  the  aftbody  upper  surface  force  in  the  x  and  z  directions,  respectively, 
and  Xeff^Zeff  are  the  components  of  the  forebody  upper  surface  force  in  the  x  and  z  directions.  For  the 
lower  surfaces,  integration  of  Eqs.  36  and  37  produce 


Fed  =  -P2U/secrpL  sin +  costal ^ 


—  Xed^  T  Z Qdk 


F gh  —  PsPnk  —  Xgtii  T  Z gfik 


(40) 


where  X^  and  Z^  are  the  components  of  the  lower  surface  forces.  Equations  39  and  40  give  the  rigid  body 
steady  forces  on  this  vehicle  due  to  pressure  distributions  on  the  upper  and  lower  surfaces  (cd  and  gh).  It 
must  still  be  determined  where  these  forces  act  on  the  surfaces.  To  do  this,  note  that  the  steady  pressures 
on  the  surfaces  are  simply  the  pressure  behind  the  shock  wave  or  expansion  fan  or  the  freestream  pressure, 
depending  on  the  angle  of  attack  of  the  vehicle. 
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V.  Afterbody 


The  flow  on  the  afterbody  of  the  vehicle  is  bounded  by  the  vehicle  surface  (surface  ef  in  Fig.  1)  and  a 
shear  layer  between  the  freestream  atmosphere  and  the  exhaust  gas  of  the  engine.  Hence,  external  nozzle 
analysis  must  be  performed  to  determine  the  pressures  on  surface  ef.  According  to  Chavez,^  the  pressure 
distribution  along  the  nozzle  surface  can  be  approximated  by 

Pefis)  «  - ^ 

1  +  - ^ 

cos(ti  {7  +  T2 

where  Pg/  is  the  pressure  on  the  afterbody,  Pg  is  the  pressure  at  the  engine  exit,  Poo  is  the  freestream 
pressure,  and  s  is  the  distance  from  the  lower  apex  (point  e)  to  the  point  of  interest  along  the  vehicle’s 
afterbody  surface  (see  Fig.  1).  The  force  produced  by  the  external  nozzle  can  be  calculated  by  integrating 
Eq.  41  over  the  rear  ramp  of  the  vehicle: 


"ef 


-I: 


=  (’'l,C/  +  T2) 


Pe 


1  + 


s 

_ l±a _ 

cos  (ti  {7+1-2) 


Performing  the  integration  and  simplifying  yields 


LaPePcx 


(42) 


(43) 


Equation  43  provides  the  magnitude  of  the  force  due  to  the  external  nozzle.  Its  direction  is  perpendicular 
to  the  rear  ramp.  Hence,  the  vector  force  due  to  the  external  nozzle  is 


Fe/  = 


_ LgPePoo 

cos  fri  rr  +  To)  (P^ 


-InA 


sin  (rin  +  T2)  i  -  cos  (ri , 


where  Xef  and  Zef  are  the  axial  and  normal  force  components  of  the  external  nozzle  force. 

Eor  use  in  calculating  stability  derivatives,  it  is  necessary  to  determine  the  force  on  the  rear  ramp  due  to 
perturbations  in  the  velocities  u  and  w  and  the  rate  q.  The  differential  force  on  the  rear  ramp  is 


dFgj  —  (  Pgj  pefaef  \y ef  '  ^ef])  dA^fW^f  (45) 

where  Pg/  is  the  afterbody  pressure  distribution  given  by  Eq.  41,  pg/,  Ug/  are  the  density  and  speed  of  sound 
on  the  rear  ramp,  Vg/  is  the  flow  velocity  on  lower  surface  ef,  ng/  is  the  normal  vector  to  lower  surface  ef, 
and  dAef  is  the  surface  element  for  surface  ef.  The  position  vector,  normal  vector,  and  surface  element  for 
this  surface  are  given  by: 


Tef  =  xi  [tan(ri^t/  +  T2)  {x  +  Xa)  —  Ptanri^t/]  k  =  xi  r^f^  k 
-Xa  <  X  <  La  -  Xa 


(46) 


ng/  =  -  sin  (ryt/  +  r2)  i  +  cos  (ryt/  +  T2)  k  (47) 

dAef  =  sec  {ti^u  T2)  dx  (48) 

Using  Eqs.  46-  48  in  Eq.  45  yields 

dFef  =  (-Pe/  -  Pefdef  {“  {u  +  grg/J  siu  (ry f/  +  T2)  +  (re  -  qx)  cos  (ryt/  +  T2)})  UefdAef  (49) 

Integration  of  the  steady  component  of  Eq.  49  yields  the  same  result  as  in  Eq.  44.  The  unsteady  components 
will  be  utilized  in  the  computation  of  stability  derivatives. 
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VI.  Control  Surfaces 


The  control  surface  is  an  elevator  located  near  the  tail  of  the  vehicle  as  shown  in  Fig.  1.  The  elevator  is 
modelled  as  a  flat  plate  hinged  at  its  midpoint  so  the  entire  surface  deflects.  The  length  of  the  elevator  is 
Le  =  17 ft.  Positive  6e  is  defined  as  trailing  edge  down.  Once  again,  the  velocity  of  the  flow  on  both  sides  of 
the  elevator  must  be  determined.  Proceeding  in  a  manner  similar  to  that  which  has  already  been  done,  it  is 
found  that 

^CL  —  (1^  COS  Se  u)  i  {  —  V4  sin  ^  ^  x  (50) 

and 

Vetj  =  (V5  COS  6e  ^  u)  i  ^  {  —  V^  sin  6e  ^  w)  k  LV  X  Ye  (51) 

where  Ve^  is  the  flow  velocity  on  the  underside  of  the  elevator,  Vejj  is  the  flow  velocity  on  the  upper  surface 
of  the  elevator,  u?  =  qj^  Tq  is  a  position  vector  from  the  vehicle  c.g.  to  an  arbitrary  point  on  the  elevator, 
V4,  V5  are  fluid  velocities  (freestream,  behind  oblique  shock,  or  behind  expansion  fan). 

The  position  vector  is  found  to  be 

Ye  =  xi-  [Zcs  +  tan ^  {x  +  Xcs)]  k  ^.^2) 

—Xcs  —  ^  cos  Se  <  X  <  —Xcs  +  ^  COS  Sq 

where  Xcs  and  Zcs  are  the  x  and  z  positions  of  the  midpoint  of  the  elevator  referenced  to  the  c.g.  As  shown 
in  Fig.  1,  Xcs  =  ~30/t  and  Zcs  =  —3.5 ft.  For  this  control  surface,  outward  pointing  normal  vectors  for  both 
the  lower  and  upper  surfaces  are  needed.  These  normals  are  computed  as 


Uc^  =  —  sin  6ci  —  cos  6ek 
UcL  =  sin^e'^  +  cos  6ek 

To  evaluate  the  cross-product  in  Eq.  51,  we  use  Eq.  52  to  obtain 

(jj  X  Yc  =  —q  [zcs  +  tan^e  {x  —  Xcs)]  'i  —  qxk 

Then,  the  differential  forces  on  the  upper  and  lower  surfaces  of  the  elevator  becomes 


(53) 


(54) 


dFeu  =  [-P5  -  Pbdb  {Veu  ■  nee,}]  ne^dAg 
dFeg  =  [-P4  -  pAttA:  {V CL  ■  nei}] 


(55) 


where 


Vejj  «  ncjj 


-  {u  —  q  {zes  +  tan^e  {x  —  ^cs)})  sin^g  —  {w  —  qx)  cos^g 
Vg^  .  ng^  =  {u-  q  {zes  +  tan  5c  {x  -  Xg^)})  sin  Sc^  {w  -  qx)  cos  5c 


and 


dAe  =  sec5edx{l)  =  sec  5cdx 

Hence,  the  upper  and  lower  forces  on  the  elevator  can  be  computed  as 


F  — 

^  eu  ~ 


Fe.  = 


/: 

£ 


[  P5  P5(^5  O^eu  '^ec/}] 
[-P4  -  PAdA  {VeL  ■  nei,}] 


■  sinJpZ  —  cos  Spk 


sec  Spdx 


sini^ei  +  cos  dgk 


sec  5gdx 


(56) 

(57) 

(58) 

(59) 

(60) 


Eq.  56  would  be  used  in  Eq.  59,  while  Eq.  57  would  be  used  in  Eq.  60.  The  steady  forces  on  the  elevator 
become 


/ 


cs  o  cos 


sin  5 ci  —  cos  5 ck 


sec  5cdx 


(61) 
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Xcs~\~  o  COS  (5e 


P4  sin  5ei  +  cos  5ek  sec  d^dx 


Evaluating  these  integrals  yields 


Yeu  =  -PbLe  -sin^i  -  cos  6ek\  =  Xe^i  +  Zeuk 


Fe^  =  —P^Le  sin^e'<^  +  COS  Sek  =  Xej^i  +  Zej^  k 


VII.  Flow  Analysis 

In  the  preceding  analysis,  the  properties  of  the  flow  on  the  upper,  lower,  and  control  surfaces  has  been 
left  general.  In  this  section,  the  properties  of  the  flow  will  be  determined.  Specifically,  the  angles  of  attack 
at  which  shock  waves  or  expansion  fans  are  created  will  be  delineated. 

By  examination  of  Fig.  1,  the  following  relationships  can  be  determined: 

if  a  =  Ti^u  freestream  :  Vi  =  Voo,  Pi  =  Poo.cii  =  ctoo 

if  a  >  Ti^u  expansionfan  (65) 

if  a  <  Ti^u  shock  {compression) 

The  wedge  angles  associated  with  the  upper  surface,  for  calculation  of  flow  properties  behind  the  shock  or 
expansion  fan,  are  as  follows: 

if  a  >  Ti^u =  a  -  Ti^u 
if  a  <  Ti^u  ^  ^u^hock  =  ““  +  ''“i.c/ 

The  above  information  is  used  to  determine  the  flow  properties  for  the  upper  surface,  namely,  Vi,pi,  and 
ai.  For  lower  surface  cd,  the  relationships  become 


if  OL  =  -TyL  ^  freestream  :  F2  =  Ko,  P2  =  Poo,  <^2  = 
if  a  >  shock{compression) 

if  a  <  expansionfan 


Lower  surface  cd  wedge  angles  are  as  follows: 


if  a  >  -ri,L  ^ 
if  a  <  -Tl  L  ^ 


a  +  ri,L 
=  -a  -  Til 


Physically,  for  the  scramjet  engine  to  work  properly,  an  oblique  shock  must  form  on  the  underside  of  the 
vehicle  to  increase  the  pressure  at  the  inlet.  This  effectively  places  a  lower  limit  on  the  angle  of  attack  and 
requires  that  the  angle  of  attack  be  such  that  an  oblique  shock  forms  on  the  underside  of  the  vehicle.  Hence, 
from  Eq.  67,  the  angle  of  attack  must  satisfy 

a  >  -TqL  (69) 

In  other  words,  is  an  absolute  lower  limit  for  angle  of  attack,  with  the  engine  on;  however,  the  engine 

will  cease  to  function  before  this  lower  limit  is  reached.  With  this  limit  in  place,  a  bow  shock  will  form  on 
the  underside  of  the  vehicle.  For  lower  surface  gh,  it  is  necessary  to  calculate  the  angle  at  which  the  shock 
exactly  impinges  on  the  point  g  of  the  engine  nacelle.  This  angle,  denoted  by  r^owshocki  is 

tan ri,L  +  /iA 


Tbowshock  =  a  +  tan  - J  (70) 

Let  the  bow  shock  angle  be  denoted  by  (3.  If  /3  >  r^owshock^  then  the  shock  misses  the  point  g  and  the  flow 
properties  on  lower  surface  gh  are  computed  using  the  flow  properties  behind  the  oblique  shock  (F2,P2,<^2) 
as  the  initial  conditions.  Then,  an  expansion  fan  forms  at  point  g.  If  (3  <  riyowshock^  then  either  the  shock  is 
on  the  lip  (point  g)  or  the  shock  is  inside  the  engine  inlet.  In  either  case,  freestream  properties  are  used  to 
compute  V3,p3,a3.  The  following  steps  are  used  to  determine  the  flow  conditions  over  surface  gh. 
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1.  Calculate  shock  angle,  /5,  from  surface  cd. 

2.  li  (3  >  Tbow shocks  an  expansion  fan  forms  at  point  g  and  the  flow  properties  behind  the  oblique  shock 
wave  are  used  as  input  to  the  expansion  fan  flow  equations  with  a  wedge  angle  of 

3.  If  (3  <  Tbowshocki  an  expansion  fan  forms  at  point  g  and  the  freestream  flow  properties  are  used  as  input 
to  the  expansion  fan  flow  equations.  The  wedge  angle  in  this  case  is  only  a  function  of  angle  of  attack. 
If  Q/  =  0,  then  the  flow  properties  in  region  3  are  freestream,  that  is,  V3  =  Voo,  Ps  =  Poo,  and  as  =  aoo- 
If  a  >  0,  a  shock  forms  at  point  g  with  the  shock  angle  computed  using  a  wedge  angle  of  a.  If  a  <  0, 
an  expansion  fan  forms  on  the  underside  of  the  engine  and  flow  properties  in  region  3  are  calculated 
using  a  wedge  angle  of  —a. 

For  the  control  surface,  the  flow  behind  the  leading  edge  of  the  elevator  is  determined  by  the  elevator 
deflection  angle  and  the  angle  of  attack.  More  specifically,  if  then  both  the  top  and  bottom  of 

the  elevator  experience  the  freestream.  Therefore,  V4  =  V5  =  Voo,  P4  =  Ps  =  Poo,  and  a4  =  =  aoo-  If 
Se  >  —a,  then  an  expansion  fan  forms  on  the  top  of  the  elevator,  while  the  bottom  of  the  elevator  experiences 
compression  and  a  shock  forms.  In  either  case,  the  wedge  angle  is  a  +  ^e-  If  the  a  shock  forms  on 

top  of  elevator  and  an  expansion  fan  is  on  the  bottom  of  the  elevator.  In  this  case,  the  wedge  angle  for  the 
shock  and  expansion  fan  is  —a  —  S^. 

VIII.  Total  Forces  and  Moments 

Having  determined  the  forces  on  each  of  the  surfaces,  the  moments  about  the  c.g.  that  each  force  produces 
must  be  determined.  To  do  this,  the  location  of  each  force  on  the  vehicle  must  be  computed.  Figure  3  shows 
the  forces  acting  on  the  vehicle.  Consider  first  the  upper  surface  with  the  forebody  and  aftbody  forces  given 


mg 


Figure  3.  Forces  Acting  on  the  Hypersonic  Vehicle. 
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in  Eq.  39.  The  forebody  force  acts  at  the  point 


'Pcff  ■■  (Y,-YtanTi,t/) 

while  the  upper  surface  afterbody  force  acts  at 

Fc/o  :  (- y,  -  {xf  +  y)  tanri.tr) 


For  the  lower  surfaces, 


2  ’ 


^  gh 


Xf  -  Lf  -  ^,L/tanTi,L  +  h, 


The  force  due  to  the  external  nozzle  acts  at  a  point  given  by 


(71) 


(72) 

(73) 

(74) 


Fef  :  {xf  -  Lf  -  Ln-  Xp^f  cos(ti,[/  +  T2)  ,L/ tanri^L  -  xp^^  sin(Ti,;7  +  T2))  (75) 

where  xp^j,  is  the  x  point  of  the  center  of  mass  of  the  pressure  distribution  on  the  rear  ramp  assuming  an 
axis  system  centered  at  point  e  with  x  axis  pointing  along  line  ef  and  the  z  axis  pointing  up.  In  other  words, 
the  center  of  mass  of  the  pressure  distribution  was  computed  in  a  local  coordinate  frame  (local  to  surface 
ef).  Then,  this  distance  was  referenced  to  the  c.g.  of  the  vehicle.  In  Eq.  75,  XF^f  is  given  by 


cos  (rpt/ +  T2)  In  ^ 

The  elevator  force  acts  at 

Fg  :  (  Xqs^  ^cs) 

Equations  39,  40,  44,  63,  and  64  give  the  axial  and  normal  force  components  of  the  forces  acting  on  the 
vehicle.  Now,  the  moments  due  to  these  forces  can  be  calculated.  Positive  moment  is  defined  as  clockwise  or 
the  direction  that  tends  to  increase  angle  of  attack,  while  negative  moment  is  defined  as  counter-clockwise. 
The  moment  arms  are  given  by  Eqs.  71  -  75,  77. 


1  - 


P  —  P 

-L  e  cx 


■  In  - 


(76) 


Mcff  — 


2 


(1  -  tan^  Ti^u) 


(78) 


Mcf,  =-PiXa  tan^  Ti^u  (^/  +  y)  +  ^ly  (79) 

Merf  =  -P2Lf  tanri,L  ^  p^f  (^Xf  -  (80) 

M,h  =  PsLn  (^xf  -Lf-^^  (81) 

Mef  =  ±Xef  [L/taUTgL  -  Xp^^  slu  (ti,[/  +  T2)]  -  Zef  [xf  -  Lf  -  Ln,-  Xp^j.  COS  (rgjy  +  T2)]  (82) 

where  xp^f  is  given  in  Eq.  76  and  Xg/,  Zef  are  expressed  in  Eq.  44.  The  sign  on  the  first  component  of  Mef 
will  depend  on  the  vehicle’s  geometry.  The  rule  used  to  determine  the  sign  is  as  follows: 


+  if  L/ tanrpL  -  sin  (rpt/ +  r2)  >  0 

-  if  L/ tanrpL  -  sin  (rpt/ +  T2)  <  0 

The  moments  produced  by  the  elevator  are 

M^ejj  —  Ts-^e  sin  ^gZgg  T  P^Le  cos  5eXes 


(83) 


(84) 
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(85) 


Mej^  =  P^Lq  sin^e^cs  ~  P^Le  COS  SqXcs 

The  total  aerodynamic  forces  and  moments  on  the  vehicle  are 

X total  =  ^cff  +  ^cfa  +  ^cd  +  ^gh  +  ^ef  P  P  ^eu  (^6) 

Ztotal  =  ^cff  +  ^cfa  +  ^cd  P  Zgh  +  Zef  +  (87) 

P^total  —  ^cff  T  ^cfa  T  P^cd  T  P^gh  P  P^ef  T  P^cl  T  P^eu 


IX.  Stability  Derivatives 


The  stability  derivatives  of  interest  for  this  vehicle  configuration  are  Cza^  Cxa^  Cmcx^  Czq^  and  Cmq- 
These  computations  are  complex  and,  in  an  attempt  to  simplify  them,  the  following  notation  will  be  used 
when  necessary: 

_  {Lg-Xaf 


-^e/2  2 

T  (Lg  —Xg)^ 

-^e/s  3 

r  (Lg  —Xg) 

^e/4  —  4 

r  (Lg  —Xg)^ 

^e/2  —  5 


+ 


3 


+  ■ 


A 


Ai  =  pette  +  Pco)(ae^aco)(L„  Xa) - ^  (pe  -  Poo)  {La  -  Xa)  -  jf-  (fle  “  Ctoo)  {La  “  Xa) 


(89) 


(90) 


(91) 


^2  —  {Pe  Poo)  P  <^oo)  ^2  {Pe  Poo)  ipe  ^oo)  {^a  ^a) 

A  _  (Pe  Poo^pe  Q-cxj) 

^3  -  Ll 

/i  =  tan(Ti,[/  +  T2) 

/2  =  Ltanrgfy 
/s  =  sin  {ti^u  +  T2) 
fi  =  cos  {tix  +  T2) 

Additionally,  the  pressure  on  the  rear  ramp  of  the  vehicle  (surface  ef)  is  given  by  Eq.  41,  which  is  a  function  of 
the  distance  moved  along  the  rear  ramp.  Thus,  this  pressure  distribution  is  not  constant  along  the  surface. 
At  point  e,  the  pressure  is  given  by  the  engine  exit  properties  {Pe),  while  at  point  f  the  pressure  is  the 
freestream  pressure  (Poo)-  An  obvious  choice  for  the  temperature  distribution  on  the  rear  ramp,  Tg/,  is  a 
form  like  the  corresponding  pressure,  so  that 

Tef{s)  ^ - T  (92) 


1 


17  +  ^2) 


(t-i) 


Then,  the  speed  of  sound  and  density  can  be  calculated  from  the  definition  of  speed  of  sound  and  the  perfect 
gas  law: 


(^ef  —  \J  'jRTef 


—  P 

ef 


Pef  —  RT 


(93) 


Both  tte/  and  pef  factor  into  the  stability  derivative  calculations.  Unfortunately,  the  expressions  in  Eq.  93, 
when  used  in  the  stability  derivative  calculations  do  not  allow  determination  of  a  closed-form  solution.  In 
order  to  facilitate  a  closed-form  solution,  the  following  approximations  for  the  speed  of  sound  and  density 
on  the  rear  ramp  will  be  used  in  the  stability  derivative  calculations: 


Pef  =  {X  -  {La  -  Xa})  +  Pe 

aef  =  (^x  -  {La  -  Xa})  +  Ug 

where 

-Xa  <  X  <  La  -  Xa 

These  are  first-order  approximations,  which  capture  the  boundary  conditions. 


(94) 


(95) 
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A.  a  Derivative  of  Z-Force  Coefficient 

To  compute  the  change  in  Z-force  component  due  to  a  change  in  angle  of  attack,  the  infinitesimal  force 
expression  that  contains  the  vertical  velocity  perturbation,  w,  must  be  integrated.  Thus, 


(96) 


where  j  {d'Fu)z-w  differential  force,  on  the  upper  surface,  in  the  z  direction  due  to  w  motion.  Using 

the  appropriate  differential  force  elements  from  Eq.  25  in  Eq.  96  produces 


Qoo^ 


nXf  nXf 

S'  I  /  - 

\J—Xa  Jxf—Lf 

PXf—Lf  pLa—Xa 

j  Xf—Lf—Lr,.  J  —Xn 


(Cz), 

1 

QooS  _ 

1 

Qoo^ 

J  Xf-Lf-L, 

nXf  nXf 

/  —piaiw  cosTi^udx  —  I  p2a2W  cosTi^Ldx 

d  —Xa 


^Xf-Lf 


pLa  Xa 


—psaswdx 


—pefaefW  COS  {ti^u  +  ^2)  dx 


Performing  the  integrations  yields 


1 


(Cz)„  =  [“/5l“l  cos  ri,[/  {Xf  +  Xa)  -  P2a2  cos  -  psasLn] 

cos  (ti  u  +  T2)  r  .  ^  ^  ^7-1 

— - ;; -  [d^lLa  +  ^2^e/2  +  ^S^e/s] 


QooS 


If  ^  1  then  «  a.  Thus,  w  =  VoO;.  Therefore, 

Too  Too 


dCz  To, 


da  goo  <5 


[-piai  cosTi,c/I/  -  p2a2  costi^lT/  -  pza^Ln] 


TooCOs(Ti,j7  +T2) 
Qoo^ 


[^1-ha  +  ^2-he/2  +  ^S-he/a] 


(97) 


(98) 


(99) 


B.  a  Derivative  of  X- Force  Coefficient 

The  change  in  X-force  due  to  a  change  in  angle  of  attack  can  be  calculated  using 

(Cx)„ = A  /  = A  r  + r 

J  d —Xa  Jxf—Lf 


1 


Qoo^ 


pXf—Lf  pJxa—Xa 

/  +  /  (dFaf) 

J  X  f—L  f—Ln.  J —Xn. 


^f-^f 

La 


'Xf-Lf-L 


Substituting  in  the  differential  forces  yields 

1 


{Cx)y,  = 

1 


Qoo^ 


Qoo^ 

X  f  —L  f  —Ly 


f  r^f 

piaiw  smri^udx  —  /  p2a2W  sinri^Ldx 

:a  dxf-Lf 


L 

—p^a^w^dx  +  J 


La  Xa 


Pefttef-W  sin  (ti,[/  +  Tz)  dx 


(100) 


(101) 


Performing  the  integrations,  simplifying,  and  assuming  that  ^  ^  1  so  that  ^  ~  a  gives 

^  tToo  [piai  sinTi,[/L  -  pzaz  sinxi^Lf]  +  [^iTa  +  AzLef^  +  AzLef^]  (102) 


da  qooS 


Qoo^ 
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C.  a  Derivative  of  Pitching  Moment  Coefficient 

For  CMai  find  the  contribution  to  the  pitching  moment  due  to  a  velocity  w: 


(Cm)„  =  ^  [y  ^  m.-n,  -  J  ^ 

2  f^f  r^f  r^f~^f 

=  /  (dFcf)^_^  +  /  {dFcd)x_yj  +  / 

J —Xa  Jxf—Lf  J  Xf—Lf—. 


I  f^f  r^f 

-  —X  /  ^  {dFcf)^_^  +  {dFcd)z-t 

J —Xa  Jxf—Lf 


{d¥gh)^_^  +  /  {d¥ef)^_^ 


/  (rfFe/), 


Substituting  the  appropriate  expressions,  performing  the  integrations,  simplifying,  and  assuming  that  ^  ^ 
1  so  that  ^  a  yields 


OCm  li  p  .  2,  •  7-21 

^  802^"^  [-piaitanryt/sinryt/L  +  tan  ryi,  smryL^/J 

+  ^  '  q  ^  (^l^e/2  +  ^2^e/3  +  AsLef^)  +  {flXa  “  /2)  (^l^a  +  ^2^0/2  +  ^S^e/s)] 

+  ^  [plCil  COSTyc/  (xj  -  +  P2a2  COSTyL^/  (2x/  -  L f)  -  p^a^^Ln  {Ln  “  2x f  +  2L/)] 


f/ooCOS  (ryt/  +T2) 


[^l^e/2  +  ^2^0/3  +  ^3^0/4] 


D.  g  Derivative  of  Z-Force  Coefficient 

The  Z-force  coefficient  due  to  a  pitching  motion  is 


1 


+  - — c  /_ 

QoO>^  Jxf — Z/ — Ln 


{dFef), 


Performing  the  integrations  and  simplifying  yields 


dCz  I  L,  ^  .  r. 

— —  =  ^Piai—  COSTi,;7  -  tanTi,;7  SinTi,[/L 

dq  QooS  2  ’ 

+^1^/9202^  [costi,l  {2xf  -  Lf)  -  tanTi,L  sinn^ii/] 

+  - — -gp^CLZ-^  [  —  Ln  +  25/  —  21//] 

+  T  ~  ■^a/2/3^1  +  Lef^fiM  +  ^ - ^^^/l/3^2  -  ■Z^e/2/2/3^2 

+  ^ ^  ■C'e/3/4^2  +  f-^^e/4  +  ^“2 - flfs^s  —  -C'e/3/2/3^3  +  Lef^fi 
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E.  q  Derivative  of  Pitching  Moment  Coefficient 

The  pitching  moment  due  to  a  pitch  rate  can  be  calculated  using 


{Cm),  =  ^  ^  ^ 


1 


qooSc 

1 

Qoo^^ 


/X  f  nXf  nXf—Lf  pLa—Xa 

idFef),_,  +  /  (rfFe,),_,  +  /  {dF,H),_,  +  /  (dFe/) 

-Xa  Jxf—Lf  Jxf—Lf—Ln  'J —Xa 

/X  f  pX  f  pX  f — Lf  pLa — Xa 

{dFcf),_,  +  /  {dF^)^_^  +  /  {dF,H),_,  +  /  (dFef) 

-Xa  Jxf—Lf  J  X  f—L  f—Ln  J —Xa 


(107) 


Substituting  the  appropriate  expressions  yields,  performing  the  integrations,  and  simplifying  produces 


OCm 


dq  QooSc 


-piUi  tan  Ti^u 


cos  TI  u 

tan  Tpc/  sin  Tpc/  —  H - 

3  6 


{-x^f  +  2x1  +  ^xlxf) 


1  2 

-P2a2tan  ri,L 


QooSc 


Li)  Lf 

tan  ri,L  sin  -X  -  cos  ri,L  -X{2xf  -  Lf) 


QooSc 


[(^l-^e/s  +  MLefi  +  AsLefs)  (/i /s  +  /1/4)] 

[{AlLef2  +  A2Lefy  +  A^LefJ  {2fff3Xa  -  2/1/2/3  “  72/4)] 

+  ^2Lef2  +  AsLefs)  (/l /s^a  “  2fif2hXa  +  flfiXa  +  /I/3)]  (108) 

--Piai  i-x)  +  +  21:3)  +  (^3  ^  ^3^ 

tanTi,L  siiiti,l^2  +  cosri,^^  (35^  -3xfLf  +i/) 

u  o 


qooSc 


-P2CL2 


QoqSc 


psds 


3  3 


qooSc 


[(^l^e/3  +  ^2^e/4  +  ^S^e/s)  (/l/s  +  f^)  +  (^l^e/2  +  ^2^0/3  +  ^sLef^)  (/l/s^a  ~  /2/3)] 


X.  Engine 

For  hypersonic  flight,  the  propulsion  system  necessary  to  produce  the  required  thrust  is  either  a  rocket  or 
a  supersonic  combustion  ramjet  (scramjet).  The  advantage  of  a  scramjet  over  a  rocket  is  that  the  scramjet 
is  an  airbreathing  propulsion  system,  thus  eliminating  the  need  to  carry  the  oxidizer  onboard  the  vehicle. 
This,  in  turn,  allows  for  increased  payload.  The  scramjet  model  used  in  this  work  is  identical  to  that  used  by 
Chavez  and  Schmidt.^  Engine  inlet  conditions  are  primarily  determined  by  the  flow  behind  the  oblique  shock. 
The  scramjet  consists  of  3  sections;  a  diffuser,  a  combustor,  and  an  internal  nozzle.  The  flow  through  the 
diffuser  and  nozzle  is  assumed  to  be  isentropic,  quasi-one-dimensional,  while  the  flow  through  the  combustor 
is  assumed  to  be  quasi-one-dimensional  in  a  constant  area  duct  with  heat  addition.  The  working  fluid  in  the 
engine  is  assumed  to  be  a  perfect  gas  with  constant  specific  heats.  There  are  two  control  variables  which 
affect  the  engine:  diffuser  area  ratio,  Ad^  and  temperature  addition  in  the  combustor,  ATq.  Figure  4  shows 
the  engine  model. 

Inlet  conditions  to  the  diffuser,  which  are  inlet  conditions  to  the  engine  module,  are  determined  from  the 
flow  analysis  of  the  vehicle’s  lower  forebody  and  vehicle  geometry.  The  lower  forebody  flow  is  turned  parallel 
to  surface  cd.  To  determine  the  engine  inlet  conditions,  the  flow  must  be  turned  parallel  to  the  engine, 
with  the  turning  angle  given  by  Oblique  shock  relations  are  used  to  determine  the  engine  inlet  flow 
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Diffuser 

Combustor 

Internal  Nozzle 

Figure  4.  Scramjet  Engine. 
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properties,  with  the  flow  properties  in  region  2  being  the  input,  and  and  being  the  outputs 

of  the  flow  calculations.  The  engine  inlet  Mach  number,  pressure,  and  temperature  are  Me-^,  Pein^ 
respectively.  By  turning  the  flow  to  be  parallel  to  the  engine,  a  force  and  hence  moment  is  imparted  on  the 
vehicle.  The  force  is  given  by 

=  'lMiP2  (1  -  COS  (ti,l  +  a)) 

Pzi^iet  =  'yMiP2  sin  (ri,L  +  a) 

This  force  acts  at  point  d  in  Figure  1.  Hence,  the  force  acts  at 


^  inlet 

Thus,  the  moment  produced  is  given  by 


Xf  —  Lf  LftaHTi^L 


Minlet  —  LftwCiTi^LFxi^i^t  {Xf  Lf) 


(110) 


(111) 


Now,  the  flow  must  be  propagated  through  each  section  of  the  scramjet  to  determine  the  engine  out¬ 
put  properties.  For  the  diffuser,  continuity  is  applied  to  determine  the  Mach  number  at  the  diffuser 
exit /combustor  inlet: 


[1  +  P7  -  1)  M2]  [1  +  -  1)  m2  J  ) 

— 


M2 


M2 


(112) 


where  is  the  Mach  number  at  the  combustor  inlet  and  is  the  Mach  number  at  the  engine  inlet.  The 
pressures  and  temperatures  at  the  combustor  inlet  are  given  by 


Pr  =  Pe. 


-l  +  p7-l)M2/ 
.  l  +  p7-l)M2  . 


T 

7-1 


(113) 


-l  +  p7-l)M2/ 
.  1+P7-1)M2  _ 


(114) 


where  Pc  and  Tc  are  the  pressure  and  temperature  at  the  combustor  inlet,  and  Pg^n  pressure 

and  temperature  at  the  engine  inlet.  For  the  combustor,  the  exit  Mach  number,  temperature,  and  pressure 
are  calculated  using 


M3[l  +  p7-l)M2]  ^  M3[l  +  p7-l)M2]  ^  M2  Tp 

(7M2  +  if  (7M2  +  if  (7M2  +  if  Te 


(115) 


Pn  =  Pc 


+  1 

7M2  +  1 


/7M2  +  1  Mn 

'v7M2  + 1  M, 


2 


(116) 

(117) 


where  M„,  Pn,  and  T„  are  the  Mach  number,  pressure,  and  temperature  at  the  inlet  to  the  nozzle,  respectively 
and  To  is  the  increase  in  total  temperature  across  the  combustor  due  to  the  combustion  of  fuel.  For  the 
nozzle,  the  exit  properties  are 


[1  +  1(7_1)M,2]('-)  _  [1  +  1(^_1)M2](-') 

M|  ^  M2 


(118) 


where  Mg  is  the  Mach  number  at  the  combustor  inlet  and  M2  is  the  Mach  number  in  region  2  (engine  inlet). 
The  pressures  and  temperatures  at  the  combustor  inlet  are  given  by 


-1+  1(7-1)  M2- 
_1+  1(7-1)  M2_ 


(119) 
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-i+i(7-i)Mr 
.l  +  i(7-l)M|. 


(120) 


where  Mg,  Pg,  and  Tg  are  the  Mach  number,  pressure,  and  temperature  at  the  engine  exit,  respectively,  and 
is  the  internal  nozzle  ratio  defined  as  the  ratio  of  the  nozzle  exit  area  to  the  nozzle  inlet  area. 

The  thrust  per  unit  width  of  the  engine  module  is  given  by^ 


T  =  mo  {Ve 


Vqo)  T  {Pe 


{P2  -  Poo) 


Ad  An 


(121) 


where  T  is  the  engine  thrust,  rfio  is  the  mass  flow  through  the  engine,  14, Pg  are  the  flow  velocity  and  pressure 
at  the  engine  exit,  14,^2  are  the  flow  velocity  and  pressure  at  the  engine  inlet,  Poo  is  the  freestream  pressure, 
Ae/h  is  the  exit  area  per  unit  width.  Ad  is  the  diffuser  area  ratio,  and  is  the  nozzle  area  ratio.  The 
moment  produced  by  the  thrust  force  is 


M 


engine 


=  L/ tanrpL 


(122) 


The  mass  flow  through  the  engine  is  a  function  of  the  shock  angle.  Essentially,  in  order  to  calculate 
rho,  it  is  necessary  to  calculate  how  much  mass  flow  the  engine  captures.  Figure  5  sets  up  the  geometry  for 
this  calculation.  It  is  assumed  that  the  vehicle  is  operating  such  that  the  lower  forebody  surface  see  either 


Figure  5.  Geometry  for  Mass  Flow  Capture  Area. 


the  freestream  {a  =  or  a  shock  forms  off  the  forebody  surface.  Let  the  engine  inlet  capture  area  be 

denoted  by  Aq  and  the  spill  area  be  denoted  by  Ag.  Then,  if  freestream  conditions  occur,  Aq  =  hi.  When  a 
shock  forms,  use  Figure  5  to  compute  the  capture  area,  which  becomes 


where  jSicd  is  the  shock  angle  for  lower  surface  cd  and  0i  is  defined  as 

0,  =  tan-i  ( +  ^ 


(123) 


(124) 
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The  mass  flow  through  the  engine  becomes 


—  Poo-^o 


RTo 


"^0 


(125) 


With  the  inclusion  of  the  stability  derivatives,  the  engine  inlet  turning  force  and  moment,  and  the  thrust 
and  resulting  moment,  the  total  aerodynamic  forces  and  moments  on  the  vehicle  are 


Xtotal  —  ^cff  +  ^cfa  +  ^cd  +  Xgh  +  W/  +  +  W[/  +  QooS  a  +  +  T 


(126) 


Ztotal  —  Zcff  +  Zcfa  +  Zed  +  Zgh  +  Zef  +  Zei^  +  Zeu  +  QooS  +  FQooS  - h  (127) 


dq  2Vo. 


FItotal  —  FIcf.  +  FIcfa  +  Fled  +  Mgh  +  Mef  +  Mej^  +  Mejj  +  QooSc—- - a  +  qooSc-- - — - h  Minlet  +  Mengine  (128) 

C/OL  ^Voo 


XI.  Results 

At  this  point,  only  initial  simulation  results  are  available.  Currently,  the  model  has  been  simulated  open- 
loop,  with  a  fixed  control  surface  deflection,  to  ensure  that  the  model  is  operating  correctly.  One  point  of 
interest  that  can  be  obtained  from  this  simple  simulation  is  the  contribution  to  the  forces  and  moments  due 
to  the  inclusion  of  unsteady  effects.  Figure  6  shows  the  steady  and  unsteady  X  force.  Figure  7  shows  the 
steady  and  unsteady  Z  force,  and  Figure  8  shows  the  steady  and  unsteady  pitching  moment.  Obviously,  the 
unsteady  components  will  have  an  impact  on  the  total  forces  and  moments,  as  these  terms  are  not  negligible 
compared  to  the  steady  terms. 


Steady  and  Unsteady  X  Force  vs.  Time 


Figure  6.  Steady  and  Unsteady  X  Forces. 
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^  Steady  and  Unsteady  Z  Force  vs.  Time 


Figure  7.  Steady  and  Unsteady  Z  Forces. 
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Steady  and  Unsteady  Pitching  Moment  vs.  Time 


Figure  8.  Steady  and  Unsteady  Pitching  Moments. 
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XII.  Conclusions 


In  this  work,  piston  theory  is  used  to  develop  a  model  for  the  longitudinal  dynamics  of  a  hypersonic 
vehicle.  In  particular,  velocities  of  flow  normal  to  the  surface  of  the  vehicle  are  used  in  a  first  order  piston 
theory  framework  to  determine  the  pressures  on  the  surfaces  of  the  vehicle.  The  pressures  are  then  integrated 
over  the  body  to  determine  the  forces  acting  on  the  vehicle.  Piston  theory  is  useful  here  because  it  allows 
the  inclusion  of  the  unsteady  aerodynamic  effects,  which  are  not  captured  using  other  techniques. 
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